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We consider an open quantum many-particle system in which 
there are dissipative processes. The evolution of this system is 
described by a kinetic equation for the density matrix. From the 
equation describing a random Markov process in this system, we 
obtain an equation for the single-particle statistical operator. This 
equation describes the evolution of a system of identical particles in 
a mean-field approximation. The equation for interacting particles in 
thermodynamic equilibrium was obtained. The distribution function 
of a system of interacting electrons in metals has multivalence 
in a certain region of wave vectors. Among many solutions one 
is isotropic. Other solutions have the anisotropy of the electron 
distribution over the wave vectors. The anisotropy arises as a result 
of repulsion and attraction between electrons. 
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1. Introduction 

The most general statistical description of any system in quantum mechanics is performed by using the 
statistical operator [1]. If the system state changes with time, then the statistical operator is a function 
of time: 

g^g{t). (1.1) 

When the system takes a random Markov process, this function can be found from the kinetic equation 
of the form [2 - 14] 

i Ti g ~ , g'^ + i Ti d g , (1-2) 

which is called the Liouville - von Neumann equation. In this equation, H is the Hamiltonian of the 
system; d is the so-called dissipative operator describing dissipative processes in the system. The statistical 
operator should be normalized, self-adjoint and positive. Any solution of (1.2) that has a physical 
significance and has these properties at the initial moment of time to should have them for all t > to. So 
the equation (1.2) should have the form 



Ihg^ [H, g] +ih ^ Cki (^[cikg, a^] + [dk, g d^] 



This equation is called the Lindblad equation [6]. Here ak is an operator, fc = 1, 2, , 
operator, C'ki - constants. 

Let us write the equation (1.3) for the quantum harmonic oscillator: 



(1.3) 

the number of the 



'.hg= [H, g] + — ihA(^[dg, a'^]-t-[a, g d^]^ + — ih B (^[d~^ g, d] + [d'^ , g d] 



(1.4) 
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where 

H = huj^a+ , (1.5) 

A and B are constants. The operator a has the form 

1 / P 



+ ^/^X . (1.6) 



Equation (1.4) surprisingly accurately describes the change of the state of the quantum harmonic oscillator 
with time, provided that 



^=e^^-, (1.7) 



where (3 = 1/kT is the inverse temperature. 



2. The hierarchy of statistical operators 

Consider the system consisting of N identical particles. The statistical operator g describing the state 
of this system can be written as follows: 

Q=Qil,2,...,N). (2.1) 

Here the numbers in parentheses denote the indices of the variables which are affected by this operator. 
The statistical operator like any other operator describing the state of the system of identical particles 
must be symmetric: 

e(..., i,..., j, ...) = g{..., j, ?,...) . (2.2) 
Let us assume the following normalization condition for the statistical operator: 

Tri .^^(l,...,Ar) = An. (2.3) 

If the system state changes over time, the statistical operator is a function of time: t: g = g{t) . When 
we can be sure that the system takes a random Markov process, this function can be found from the 
equation 

in^ = [n, q] +indQ. (2.4) 

The Hamiltonian operator in this equation can be written as the sum 

N N N 

n{l, N) = ^^(z) + y E E ^(^' J') ' (2-5) 

i=l i=l j=l 

where H{i) is a single-particle Hamiltonian, i.e. the energy operator of a single particle without considering 
its interaction with other particles; H{i, j) is the operator of the interaction of two particles. A single- 
particle Hamiltonian may contain dissipative terms. A two-particle Hamiltonian should be symmetric: 

H{i,j) = H{j,i). (2.6) 

Under this condition, the Hamiltonian (2.5) will also be symmetric. 

Suppose that the action of the dissipative operator d on the statistical operator g leads to the expression 

dg= ^ (^[ag, d+] + [d, ga+]y (2.7) 
Here the indices k and / near the operator d are omitted . 
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Suppose that the dissipative processes occurring in the system are caused by the stochastic interaction 

of particles with a heat reservoir. Moreover, each particle interacts with it independently of other particles. 
In this case the many- particle operator a in the expression (2.7) can be written as the sum 

N 

a{l,...,N)=Y,a{i), (2.8) 

i=l 

where the operator a{i) characterizes the effect of the heat reservoir on one of the particles. 
Substitution of the sum (2.8) in the formula (2.7) gives 

N N 

^^=2^ E( S' hh] +[a{i), Q a + (i)] ) . (2.9) 

i=i j=i 

It is impossible to solve the equation (2.4) with such a complex right-hand side. Therefore the practical 
interest is carried by the single-particle kinetic equation for the statistical operator ^(1), which is defined 
by the relation 

^^^^ = jN~[y. ^' • ^^-^^^ 

This definition and condition (2.3) implies the normalization condition 

Tri^(l) = Ar. (2.11) 
We define the two-particle statistical operator ^(1, 2) as 

g(l:2)^ Tr3 ,^g(l, 2,...,^). (2.12) 

This operator satisfies the normalization condition 

Tri2^(l,2) = Ar(7V-l). (2.13) 



3. The quantum kinetic equation for single-particle density matrix 

To obtain the equation for the single-particle operator q{1), let us apply the operation of coagulation 
Tr2... jv to both sides of the equation (2.4). As 

Tr, [%), ^(l,...,iV)] =0 

by definition (2.10), we obtain 

JV 

Tr2...jv E [^(0' ^(•••' ^' •••)] = (^ - 1)! [H{1), m] ■ 

i=l 

Due to the identity 

Tr,Jij(^,j), ^(1,..., A^)] ^0, 
taking into account the property (2.6) and the definition (2.12) we have 

N N 

i=i j=i 



(JV N N 

E [^(1, i), ^(1, A^)] + y E E [^(^' J')' S{h N)] 
j=2 i=2 j=2 




= {N-iy. Tr2 [^(1,2), ^(1,2)] . 

Similarly 

Tr2...jv dg = I Tr^...^ (^[ail) g, a + ( 1 )] + [ a(l ) , ^ a + ( 1 )] + 

AT 



^ d + {l)] + [a{i), ga + {l)] + [a{l) g, a + (i)] + [a(l), ga + {i)]) + 

1=2 

N N \ 
i=2 j=2 / 

l(iV -1)! ([«(!) ^(l),« + (l)] + [a(l), ^(l)a + (l)]+[Tr2a(2)^(l,2), a + (l)] 



+ [a(l), Tr2^(l,2)o + (2 



as 



Tr, [a(2), ^a + (j)]= 0, Tr^- [a{i) g, a + {j)]=0. 

Bringing together these expressions, we obtain the equation 

in§{l) = [^(l),^(l)] +Tr2 [H{1,2), g{l,2)] + 
+ 1 i ?i([a(l) ^(1), a+(l)] + [a(l), ^(1) a + (l)] + [Tv^ a(2) ^(1, 2), a+(l)] + 

+ [a(l), Tr2^(l,2)a + (2)]). (3.1) 
We denote the matrix elements ^(1), ^(1, 2), H{1), H{1, 2) and a(l) in an a-representation as follows: 

£•010; = giV ) ga^a2,a[a'^ = £'12,1'2' , 

Ha^a[ = Hiy , H^^a^^a'^a'^ = -^12, 1'2' , (la^a{ = j 

and write the equation (3.1) in matrix form as: 

i h gii> = ^^(i?i2 £'21' — ^12 ^^21' ) + ^^(^^12,34 £'34, 1'2 — ^12,34 ^^34, 1'2 ) + 



"2 as 04 



+ ^*^ 2IE(2'*12^23«31' -«l2«23e31' " ^12 43«31') + 
+ EEE(^12,34 (a42«tl' -«42«31') + ( « 12 «34 " «|2 «34 ) £-42, 31' ) • (3-2) 



a2 013 at 



The two-particle density matrix describing the fermion system state must be antisymmetric. This 
matrix can be approximately expressed by means of the single-particle density matrix as follows: 

£'l2,l'2' ~ ^11' £"22' ~ ^12' £"21' • {'^■'^) 

Substituting this expression into the equation (3.2), we obtain the equation for the single-particle density 
matrix 

iTig-lV = g2V - ^12^^21') + EXK^'^lZ ^23 «31' - «l'2«23 ^31' " ^12 a23 «31' ) + 

+ 1^^ mm(^'i2fi'34 (a43«2i' -"43021' - a23a4i' +43041') + 

CX2 Ct3 OiA 
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+ («12 "34 - «12 0-34 - tt^j + «14 «32) £"43 £'21') ' (3-4) 

where 

Hi2 =Hi2 + 2j2Yl ^13, 24 ^43 • (3.5) 



as cn 

Let us write the equation (3.4) in operator form 



iTiQ= [H, q] +^ih(^[aQ, o+] + [a, ^ a+] ) + 

+ ^in([Tr2(a+^)a-Tr2(o^)o+, ^] + [Qa+ , go] + [a+ q, a q]^ , (3.6) 
where ^ = ^(1), 

F = ff(l) + 2 Tr2 [^(1,2), ^(2)] (3.7) 

is an averaged single-particle Hamiltonian. The equation (3.6) can be given more compact form if we use 
the notation 

H =H+^ih{Tv2{a+Q)a-Tv2{aQ)a+^ . (3.8) 

In this case we have 

int=[H',Q] +^ih[[aQ, a+{l-Q)] + [{l-Q)a,Qa+]) . (3.9) 

In general the right sides of equations (3.8) and (3.9) may contain several terms describing dissipative 
effects which depend on different operators dk- Taking this into account we generalize the equation (3.8) 
and write it as follows: 

inQ=[H , q] +in^Cki{^[ahQ,al{l- q)] + [{1- Q)ak, Qcif]^ , (3.10) 

k,l 

where 

H =H+ in ^ Cki (Tr2(a+^)afc-Tr2(afc^)a+) . (3.11) 

k,l 

We write the equations (3.10) and (3.11) in matrix form: 

ihQiv = X!(-f^l2 521' - ei2i^21') + ^12,31' £'23 -r23,12 P3I' " r3i/,23£'l2) + 

+ ^ ^ ^ 511 X](^12 £'34 (^43^ 21' ~ ■'^21', 43 ~ ■'^23,41' + ^41/^ 23) + 
0:2 as 04 

+ (^12,34 ^ ^34^2 ~ r;^4_32 + 14) £143 ^21' ) ' (3-12) 

= ifl2 + ^ ^ (2 ifi3,24 + \in (ri2,43 ' ^43,12)) ^43 , (3.13) 

013 ai 

ri2,3i' =2 X]<^fc'"'i2,fea+, (3.14) 



where 



4. The principle of detailed balance 

Suppose that at some point in time the density matrix becomes diagonal: 

Qw = Wi 6n' , 



(4.1) 
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where Wi = Wc^ is the probabihty of the state ai being filled by a particle. The substitution of matrix 
(4.1) into (3.12) leads to the equation 

Wi = Y^ {Pi2 (1 - VFi) - -P21 (1 - W2) WC) , (4.2) 

where 

^12 = r 12, 21 (4.3) 

is the probability of the transition of a particle from state 0:2 to state ai during a time unit. The transition 
probability P12 can always be represented as 

Pi2 = P[^ e-i^^'''-^'\ (4.4) 

where 

p(o) _ p(o) 

-'^12 — -'^21 I 

Si - the average energy of a particle in the state at, which is an eigenvalue of operator (3.13): 

£1 =£1+2 ^i?12,12 W2, (4.5) 

where Si — Hn. 

When a fermion system comes to the state of thermodynamic equilibrium, the particle distribution 
over the states must be subjected to the principle of detailed balance. We show that this distribution 
follows from the equation (4.2). Since the probability Wi no longer depends on time, the right side of 
this equation is equal to zero. Moreover, according to the principle of detailed balance each term should 
be equal to zero. Thus the principle of detailed balance in this case is expressed by 

Pl2{l-Wi)W2= P2i{l- W2) Wi . (4.6) 
We substitute the expression (4.4) into this equation. After simple transformations we obtain the equality 



Wi W2 

in which the left side depends on ai, and the right one - on a2- This is only possible if both sides are 
equal to the same constant value. Denote this value e~^^, where fj, is the chemical potential. We get the 
equation 

l^ = ef^(^-^). (4.7) 

Using the equality (4.5) we could show that this equation has an anisotropic solution [15]. 

From the equation (4.7) we find that the equilibrium distribution of noninteracting particles over states 
is described by the Fermi - Dirac function 

W = ^7 ^ . (4.8) 



5. The electrons in a metal 

Consider electrons in metal. Equation (4.7) is transformed to 

In = ^ (ej^ _ ^) . (5.1) 

Here, the function describes the electron distribution over wave vectors k, 

£k = £k + Yl ^kk' (5.2) 
k' 
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is the average energy of an electron, is the energy of an electron without taking into account its 

interaction with other electrons, ej^j^' is the interaction energy of electrons with wave vectors k and k'. 
Assume the following approximate formula for the interaction energy 

^kk' = -'^'^k+k' - "^^k-k' I (5-3) 

where I is the repulsion energy of electrons with vectors k and — k', J is the attraction energy of the 
electrons with vectors k and k'. The substitution of (5.3) into (5.2) gives 

Sk = £k + Iw-k- Jwis., (5.4) 

and the equation (5.1) becomes 

In i — = ^ (gj^ + _ _ ^) . (5.5) 

Replace in this equation the vector k by the vector — k: 

In— — = (g _[_ j^; _ J _ „^ (5.5) 

where 



E-k = Ek • 

The equations (5.5) and (5.6) form a system with two unknowns Wy^ and w_k- 

This system has the solutions of two types. One of them describes the isotropic distribution of electrons 
over wave vectors and the other describes the anisotropic distribution. If w_k = ^k, then each of the 
equations (5.5) and (5.6) turns into the equation 



Ini^ 



= ^(^e^ + {I-J)wi,-i,) . (5.7) 



In the equations (5.5) and (5.6) the unknown functions and w_k are complex functions, in which 

the role of an intermediate variable is played by the kinetic energy of an electron: w;_k = ''-^'li^k) and 
Wk = w'2(£k)- The functions wi = wi{s) and 102 = W2{£) are the solutions of the equations 

In = — 

Wl T 



, 1-W2 

In 

where 



2e+{l- f)w2-{l + f)wi) , (5.8) 
= ^(^2e+{l-f)wi-{l + f)w2), (5.9) 



e-ii 46 
e = -r—r , r = 



J+I ' J+/ ' 

the ratio of the energy / and J is determined by the parameter 

J- I 



f = 



J + I 

Without loss of generality, we may accept the condition 

w^i(e) < W2(e). (5.10) 
Let us introduce new variables d and s by means of relations 

W2 — Wl = d , lUi+W2 = l + S. (5-11) 
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Due to (5.10) the quantity d is non-negative: d> 0. The maximal value of d is equal to one: e [0, 1]. 
The quantity s takes on values from — 1 to 1: s e [ — 1, 1]. Let us solve equations (5.11) with respect to 
wi and W2- 

wi = |(l + s-d), W2 = ^{l + s + d). (5.12) 

We transform equations (5.8) and (5.9) using (5.12) to a form convenient for their numerical solution. 
To do this we first subtract one equation from another, then we add these equations. As a result, we 
obtain the following system: 



(l + rf)'-s^ 

(l-d)2-s2 



T . (l-.s)2-d2 1 



^=8 (l + :s)2-rf2 + 2^^ + ^)^- 



(5.13) 



These equations allow us to represent the energy e and the probabilities Wi and ■u;2 as functions of the 
parameter d. Using these functions it is not difficult to construct the graphs of the functions wi = wi(e) 
and W2 = W2{e) for different values of temperature. These graphs arc shown in Fig. 1 for the case when 
the interaction energy of electrons / and J is such that J = 2)1. At the same time f = 1/2. 

The solution of equation (5.7) is also a function of electron kinetic energy £\^: t«k = WQ{e\^). The 
function wq = Wo(e) can be found from the equation 

ln^^ = l(e-/«;o), (5.14) 

Wo T ' 

which is a corollary of equation (5.7). The solution wttq = w^o(e) of equation (5.14) is among the solutions 
of equations (5.8) and (5.9). Indeed, if we put in these equations wi = W2, then each of them takes the 
form (5.14). 

Fig. 1 shows the graphs of all three functions wq = WQ{e), wi = wi(e) and W2 = W2{e). From this 
figure it is clear that at sufficiently low temperatures the distribution function for some values of the 
energy e can take not one but several values. The multivaluedness of functions w = w{€) suggests that at 
the same temperature there can be different equilibrium macrostatcs of a system of conduction electrons 
in metal. These macrostates differ from each other by the distributions of electrons in Bloch states. In 
fact, only the macrostate of electrons, in which their energy is minimal, is realized, provided that this 
state is stable and cannot be destroyed by any external effects. 

The graphs in Fig. 1 give an idea of the changes in the distribution of electrons in states that occur 
when the metal temperature changes. When the temperature r > 1, only the isotropic distribution of 
electrons over wave vectors described by the function ru]^ = wo(£k) is possible. The graph of the function 
Wo = Wo{s) for all values of the temperature passes though the point f2 with coordinates e = /x + ^ ( J — /) 
and w = |. As the temperature decreases, the slope of the curve at this point increases. When the 
temperature is sufficiently low, the curve of = wo{e) bends so that it becomes similar to the letter Z. 

For T = 1 on the curve wo = wo{e) at the point a closed curve Z originates. Its size increases as the 
temperature decreases. The shape of the curve also changes. The value r = 1 corresponds to the critical 
temperature 

r. = 4±f. (5.15) 

At T € (r', 1), where r' is some critical value, the vertical line intersects the curve Z no more than in two 
points (the curve 2 in Fig. Ic). At r < r' the curve Z bends so that the vertical line intersects it four times 
in some places (Fig. lb). In the limit, r ^ 0, the curve Z transforms into a polygon AB1C1OC2B2A, 
that looks like the letter Z. This polygon is shown in Fig. la. The curve Wo = Wo{s) divides the curve Z 
into two parts. The lower part corresponds to wi = Wi{e) and the upper - to W2 = W2{e). 

Let us give more detailed consideration to the distribution of electrons over states at T = 0. In the 



- 9 - 



limit, r — )• 0, the solution of (5.14) takes on the form 

1 if e<n + J- I , 
wois) = { "7 — ^ if n<s<iJ, + J — I, 



J -I 

\{ £>ll. 



(5.16) 



The graph of this function is shown in Fig. la as the dotted broken line B2AOC1. 

In the limit, r — ^ 0, the equations (5.8) and (5.9) give the following dependence of the probability w on 
the electron kinetic energy £, which describes the anisotropic distribution of electrons over wave vectors: 



w{e) 



(5.17) 



1 if £<pL-I, 

Wile) if fJ- — I<£<IJ, + J, 
if s> n + J , 

where i = 1 or 2. The values of the functions Wi = iui{e) and W2 = W2{s) form pairs, such that 

wi{e) = and W2{e) = 1 (5.18) 

at /i — 7 < e < + J, or 

iui{s) = — {e — fj, + 1) and W2(e) = 1 (5.19) 

at/i — / <£</i+J — 7, or 

«;i(£)=0 and W2{s) = ^{e - jj,) (5.20) 
at < £ < /i + J. The broken line ABiCiO in Fig. la corresponds to wi = w-i{e), and the line AB2C2O 

-to W2 = W2{e). 



6. Conclusions 

The quantum Markov equation for iV-particle density matrix describes the evolution of a strictly open 
system with allowance for the dissipative operator. From this equation, we obtained the approximate 
kinetic equation for a single-particle density matrix. This equation describes the evolution of a system of 
identical particles. Wc obtained the integral equation for the function of distribution of particles in states 
where a many-particle system is in thermodynamic equilibrium. We showed that for equilibrium systems 
of interacting particles the distribution over wave vectors is anisotropic. 
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Fig. 1. The distribution function of conductivity electrons in energies for the case when 

J=3I, and at different temperatures t: a) t = 0; b) t = 0.5; c) 1 - t = 0.75; 
^-T = 0.95. 



